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Abstract. 

The anisotropy of magnetophydrodynamic turbulence is investigated by using solar wind 
data from the Helios 2 spacecraft. We investigate the behaviour of the complete high-order 
moment tensors of magnetic field increments and we compare the usual longitudinal structure 
functions which have both isotropic and anisotropic contributions, to the fully anisotropic 
contribution. Scaling exponents have been extracted by an interpolation scaling function. 
Unlike the usual turbulence in fluid flows, small-scale magnetic fluctuations remain anisotropic. 
We discuss the radial dependence of both anisotropy and intermittency and their relationship. 



The idea that small scales were statistically isotropic to a high degree, in spite of the 
anisotropy of large scales, is at the heart of the concept of universality underlying scaling be- 
haviour in turbulent fluid flow [1]. As experiments became increasingly sophisticated, evidence 
for the presence of anisotropy for both scalar fields [2] and velocity fields [3] has grown up. 
To what extent scaling laws extracted from experiments are approximated by the scaling laws 
obtained from (isotropic) cascade models [4], became an urgent issue. A workable benchmark 
is represented by the following conjecture [5,6]. Suppose that for a given statistical measur- 
able quantity, scaling features can be splitted into an isotropic and an anisotropic part. We 
might ask to what extent the ratio between the anisotropic and the isotropic part vanishes 
as the scale goes down. Of course if we are able to define quantities for which the isotropic 
part is exactly zero, we can measure scaling laws related entirely to the anisotropic part. A 
formal analysis to reveal the isotropic and anisotropic contributions to structure functions 
can be performed by considering their SO (3) decomposition [7]. This is based in a foliation of 
various j-sectors of the linear space generated by the structure functions. For a recent review 
on anisotropy in turbulent flows see [8]. 
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Anisotropy of turbulence in charged fluids, described by Magnetohydrodynamic (MHD) 
[9-11], can be due to a background magnetic field B , or a magnetic field at large scales. Since 
B cannot be cancelled out by a Galilean transformation, its direction becomes privileged with 
respect to other directions. Then MHD turbulence is anisotropic to a very large extent. A 
simple shell model for anisotropic MHD [10] showed that even if the energy injection rate is 
isotropic, or isotropy is restored at large scales, intermediate scales and mainly small scales 
are strongly anisotropic and the return-to-isotropy at small scales does not hold in MHD. This 
important result has been confirmed both by numerical simulations [11], and more recently 
by analysis of Ulysses spacecraft in fast solar wind [12], where it is shown that strongly 
intermittent and anisotropic events persist even at high frequencies. 

The study of the nature of the low-frequency magnetic fluctuations in the interplanetary 
space, in the frequency range 10 -7 Hz < / < 1 Hz, is important because it represents the 
only way to understand the behaviour of MHD turbulence in a frequency domain which is not 
accessible to any laboratory on Earth [13,14]. Anisotropy of MHD fluctuations in the solar 
wind, due to the background magnetic field transported away from the sun by the plasma 
flow [15, 16], has been investigated through the minimum variance analysis [17]. This method 
gives informations about the energy distribution among the different components of magnetic 
fluctuations. It is based on the determination of eigenvalues and eigenvectors of the variance 
matrix C Q( 3 = (B a Bp) (brackets mean ensemble averages). During the so-called Alfvenic 
periods [14] one of eigenvalues is much smaller than the other two A3 -C A2 < Ai, for example 
the ratios are Ai : A 2 : A 3 = 10 : 3 : 1 [16]. The minimum variance direction turns out to be 
almost parallel to B , so the magnetic fluctuations lie in a plane approximately perpendicular 
to B . These old studies have been used to investigate the spectral properties of the second- 
order tensor components, and to investigate the radial evolution of anisotropy. Here we 
investigate the scaling properties of anisotropy of the interplanetary magnetic field, by using 
the data analysis technique described extensively in Ref.s [5,6] and [3]. This technique allows 
us to define the scaling properties of anisotropy of the interplanetary space fluctuations, as far 
as scaling exponents of high-order tensors are considered, so that intermittency is considered. 

The data we use in this work are 6 second measurements of the Interplanetary Magnetic 
Field (IMF), as recorded by Helios 2 flux-gate magnetometer. Velocity measurements from 
Helios 2 instruments are only available at lower resolution so that the statistical accuracy 
does not allow to perform the analysis presented in this work. The spacecraft followed an 
ecliptic orbit with perihelion and aphelion located at about 0.3 and 1 AU, respectively. The 
peripheral speed of the spacecraft along its orbit was much less than the solar wind speed so 
that the latter can be still be considered aligned to the radial direction from the sun. Helios 
2 magnetic field measurements are given in a spacecraft centered reference system whose x 
axis is along the radial direction and positive towards the sun, the y axis lies on the ecliptic 
plane and is positive in the opposite direction with respect to the orbital motion, and the z 
axis is perpendicular to the ecliptic and completes the right-handed reference system. In this 
reference system, field differences at different lags r along the x direction can be obtained 
as AB a (r) = B a (x + r) — B a {x). The period of time used in this analysis spans across the 
first 4 months of 1976 when the spacecraft moved from 1 to 0.3 AU. During this period, for 
which the solar activity was near the minimum, the solar wind was clearly characterized by 
alternating fast streams and slow wind (cfr. e.g. Ref.s [13,14]). The different turbulent state 
of the fluctuations observed within these samples of fast and slow wind and the strong radial 
evolution experienced by these fluctuations (cfr. e.g. Ref.s [13,14]) convinced us about the 
necessity to analyze these data separately per wind speed and radial distance [18], avoiding 
also to include velocity shear interfaces and the only shock that was present in this particular 
data set. As a matter of fact, the radial evolution is crucial for anisotropy studies, both 
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because of the natural evolution of intermittency [19,20], and because of the rotation of the 
mean field along the Parker spiral. Thus, two different samples for each wind type have been 
selected, at two well separated distances from the Sun, namely at 0.9 AU and 0.3 AU (1 AU 
~ 1.5 x 10 8 km), each data sample including about three days of records (~ 4.3 x 10 4 data 
points). The switch between space and time sampling is obtained from the Taylor hypothesis, 
so that x = —t(v x )i the average speed being computed for each stream. 

In order to measure the contribution from the anisotropic part of the field differences, the 
whole independent elements of the n-th order tensor, computed from time series, must be 
studied 

S ai , a2 ,-,a n {r) = (i B ai{x +r) - B ai (x)] [B a2 (x +r) - B a2 (x)] ■ ■ ■ [B an (x + r) - B an (x)}) 

(1) 

For example, when n = 2 and ai = «2 = i the usual longitudinal second-order structure 
function is recovered. The latter have been studied in the past, together with the longitudinal 
higher order moments [14, 19]. For each scale, some of the components of the tensor can be 
proven to have contribution from the anisotropic part of the field only (we will use in this work 
S xz for the second order, S xxxz and S xzzz for the fourth order structure functions) [5,6,12]. 
It is thus possible to compare the scaling behavior of such components with those of the 
longitudinal structure functions (namely S xx and S xxxx ), which includes both isotropic and 
anisotropic contributions [5,6,12]. 

As a first analysis, let us discuss the scaling of the normalized longitudinal, transverse, 
and off-diagonal components of the tensor mentioned above. Figure ^ displays the scaling of 
the second order longitudinal, transverse, and off-diagonal components, for both fast and slow 
wind, at 0.3AU and 0.9AU. Each variable has been normalized to its own standard deviation 
before computing the structure functions. 

Even at a first glance, it seems evident that the diagonal elements are more or less of the 
same order at all scales. In fact, it is possible to note a scaling from the large scale values 
(Saa = 2 for uncorrelated fields) toward zero, for almost all datasets, with at most a ratio 
of order 3 for all scales. Those differences are more evident far from the sun for fast wind, 
while for slow wind they are already present at 0.3AU. The off-diagonal elements are about 
one order of magnitude smaller than the diagonal elements. For isotropic fields, they should 
vanish at large scale because of decorrelation. There are indications that this is indeed the 
case, as can be seen looking at the rightmost points of the plots in Figure ^ an< i taking into 
account the "large scale" values reported in Table [l] (see below). However, most of the mixed 
components show the same scaling as the diagonal elements, that is they do not seem to 
vanish at small scales faster than the diagonal elements. In particular, in most of the cases 
the scaling at small scales is more or less the same for off-diagonal and for diagonal elements, 
so that the anisotropic contribution could eventually survive at small scale. This can be seen, 
for example, by looking at fast wind close to the earth, where S xy decreases as slowly as the 
diagonal S aa — 2, while S xz decreases more rapidly. This important result indicates that 
isotropy, in fact, is not restored at small scales. In slow wind, the off-diagonal elements seem 
to be more scattered and, in general, slightly higher than in fast wind. This would mean that 
fast wind remains more anisotropic than slow wind. Concerning radial evolution, it might 
be said that at short distance from the sun the anisotropic contributions (off-diagonal terms) 
are less important, while they increase approaching the Earth. These observations suggest 
that perhaps the anisotropy may be hidden by the Alfvenic nature of turbulence, which is 
more evident in fast wind and near the sun (see for example Ref.s [14, 19]), while the angle 
between large scale magnetic field (the Parker spiral profile) and mean flow is of course more 
relevant as it is increased, far from the sun. This behaviour is even more evident by looking 
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Fig. 1 - The second order structure function S a p(r) of interplanetary magnetic field within fast (top 
panels) and slow (bottom panels) wind, both at 0.9 AU (left panels) and 0.3 AU (right panels), 
plotted versus the scale in log-log axes. The structure functions were computed using for each scale 
the fluctuations of magnetic field normalised to their standard deviations, so that for a large scale 
uncorrelated field the value of second order structure functions should be 2. 



at the fourth order moments, while for the odd orders moments the scaling is unclear and no 
information can be inferred (not shown) [21]. 

At variance with Ref. [12], since S ail a 2> ... ian (r) display multiple power-law scaling, we 
prefer to fit the structure functions by the following Batchelor relation [5,6,22] 



A ai ,.., an r) n {r/vT 



l + B au ... tan (r/rtf 



c„ 



[1 + D au ... <an (r/Lo)] 



2C a 



(2) 



The free parameters in equation J2J include the main scales of the flow, which need to be 
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Table I - The typical scales of the flow Lq and rj (see text), for fast (F) and slow (S) wind, at 
R = 0.3AU (©0.3) and R = 0.9AU (@0.9). 



Stream 




L (10 6 Km) 


rj (Km) 


Fast Wind R = 


0.3 


5.0 


578 


Fast Wind 7? = 


0.9 


11.9 


1920 


Slow Wind R = 


0.3 


14 


150 


Slow Wind R = 


0.9 


16.3 


2520 



estimated independently. In particular, rj is the dissipation scale, which can be estimated, 
for each stream, as the IMF typical cyclotron scale, and computed directly from the data 
as X] = 2irm p Vth/e-BQ, m p being the proton mass, V t h the thermal speed and e the electron 
charge. Note that this scale represents in fact lower bound for the actual dissipation scale, 
which cannot be directly measured from the data because of resolution limits. The integral 
scale Lq is also obtained, for each stream, from the isotropic magnetic power spectra E(k) 
as Lq = J E(k)k~ 1 dk/ J E(k)dk. Values of Lq and rj are reported in Table |U The integral 
scale Lq increases with the distance for both types of wind, and is larger for slow wind. This 
seems to indicate that turbulence is more developed in slow wind. Since the integral scale is 
growing moving away from the sun [13,14], but at the same time the estimated dissipation 
scale increases more rapidly, the inertial range and thus the Reynolds number get smaller 
moving away from the sun, as happens in decaying turbulence [4], and in disagreement with 
previous analysis [23]. 

The multiple power-law behavior of the structure functions is described by equation 
through the exponent C ai a2 ... iQri , which are the scaling parameters obtained from the fit of 
the data corresponding to the inertial range scaling, such that the usual scaling exponent can 
be recovered as £ ai) . ..,<*„ = n — 2C ai ,.... an [5,6]. This can easily be tested by dividing the 
structure functions by the factor [1 + D au ... j<Xn (r/L )] «i an d observing the power- 

law behavior, with exponent C ai ,...,a„, of the resulting extended structure function. FigureEl 
shows the behaviour of some second order mixed structure functions, for both fast and slow 
wind, and at two heliocentric distances, together with the fit with equation ©, as well as the 
compensated function and a power-law fit for one case. It is evident that equation © can 
reproduce the data with very good agreement, so that the values obtained for the reduced x 2 
are always less than unity. The most striking differences between structure functions computed 
for fast and slow wind are shown by the data samples recorded close to the sun (0.3 AU). In 
addition, a remarkable radial evolution is clearly visible, and it is stronger for fast wind, as 
already observed in different contexts [13, 14]. 

The values of the exponents obtained from the fit of the second and fourth order structure 
functions tensor longitudinal and fully anisotropic components with equation (J2J) are listed in 
table [H] for both fast and slow wind, at 0.3 AU and 0.9 AU. For both slow and fast wind 
the differences between the longitudinal and the anisotropic components are small. Once 
more this implies that the anisotropic contribution does not vanish at small scales, but its 
contribution to the scaling laws is, at least, decreasing as slowly as the isotropic part. 

In conclusion, we have investigated anisotropy of magnetic fluctuations as they are mea- 
sured in the interplanetary space, by using a method which allows to extract anisotropic 
exponents by avoiding mixing with the isotropic sector [5,6, 12]. The normalized isotropic 
longitudinal components, and the fully anisotropic components of the second and higher or- 
ders structure function tensor are shown to have similar scaling toward the small scales, 
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Fig. 2 - The second-order longitudinal structure function S X x(r) of interplanetary magnetic field 
within fast and slow wind at 0.9 and 0.3 AU, plotted versus the scale r in log-log axes, together with 
one example of the compensated function (see legend). Lines represent the fits with the function 
while the compensated function has been fitted with a power-law with exponent 0.67 ± 0.02. 

where anisotropy effects could be expected to vanish. This can be interpreted as an evidence 
that small scale anisotropy is present. Then, following Kurien and Sreenivasan [5,6] scaling 
exponents have been extracted by using a scaling function that interpolates the high-order 
tensor components. The comparison between the scaling exponents of the longitudinal and of 
the anisotropic components shows once again that small-scale fluctuations remain anisotropic. 
Thus, the return-to-isotropy invoked in the usual fluid flows does not work in MHD turbulence 
to a very large extent. 

* if if 

We are grateful to Alain Noullez for useful discussion and to an anonymous referee for 

Table II - Scaling exponents C, ai ---a n of the longitudinal and fully anisotropic components of the 
second (marked with *J and fourth order structure function components for fast (left columns) and 
slow (right columns) solar wind streams, at two different distances R from the Sun. 





FAST WIND 




SLOW WIND 




Cd 1 ■ ■ ■ OL n 


R = 0.3 AU 


R = 0.9 AU 


R = 0.3 AU 


R = 0.9 AU 




0.67 ±0.02 


0.66 ± 0.03 


0.60 ± 0.03 


0.60 ± 0.02 


Cyy 


0.69 ±0.1 


0.62 ±0.03 


0.68 ±0.03 


0.62 ±0.02 


Czz 


0.70 ±0.1 


0.66 ± 0.04 


0.68 ± 0.03 


0.57 ±0.02 




0.68 ±0.1 


0.66 ±0.03 


0.64 ± 0.03 


0.60 ± 0.02 




0.53 ±0.02 


0.67 ±0.03 


0.57 ±0.02 


0.59 ± 0.02 


(yz(*) 


0.71 ±0.09 


0.59 ±0.02 


0.67 ± 0.03 


0.62 ± 0.02 


^XXXX 


1.68 ±0.01 


1.19 ±0.03 


1.31 ±0.01 


1.18 ±0.04 


C,xzzz{,*) 


1.51 ±0.1 


1.54 ±0.03 


1.43 ±0.04 


1.33 ±0.03 


C,x raz(*) 


1.52 ±0.1 


1.55 ±0.03 


1.39 ± 0.04 


1.37 ±0.02 
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helpful comments. 
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